Nonlocal Effects and Shrinkage of the Vortex Core Radius in YNi 2 B 2 C Probed by /LtSR 
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The magnetic field distribution in the vortex state of YM2B2C has been probed by muon spin 
rotation (/iSR). The analysis based on the London model with nonlocal corrections shows that the 
vortex lattice has changed from hexagonal to square with increasing magnetic field H. At low fields 
the vortex core radius, p v (H), decreases with increasing H much steeper than what is expected from 
the \/~H behavior of the Sommerfeld constant j(H), strongly suggesting that the anomaly in 'y(H) 
primarily arises from the quasiparticle excitations outside the vortex cores. 
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The recent studies of the flux- line lattice (FLL) state in 
ordinary s-wave superconductors have revealed that the 
electronic structure of vortices is much more complicated 
than that of a simple array of rigid cylinders contain- 
ing normal electrons. One of the unexpected phenom- 
ena within this conventional model is the non-linearity 
in the magnetic field dependence of the Sommerfeld con- 
stant 7(i£) (electronic specific heat|-coefficient) observed 
in CeRu 2 EJ, NbSe 2 B, and YNi 2 B 2 C0. According to the 
above simple model where the quasiparticle excitations 
are confined within the cores of vortices (with a radius £) 
in s-wave superconductors, one would expect that j(H) 
is proportional to the number of vortices per unit cell 
and thus to the applied magnetic field H. However, ex- 
periments have revealed tifflt this is not the case for any 
of the above compoundsBB Instead, they find a field de- 
pendence like "/(H) oc y/H which is expected for d-wave 
superconductors having more extended quasiparticle ex- 
citations along nodes in the energy gap. The recent study 
on the effect of doping in YNi 2 B 2 C and NbSe 2 indicates 
that the anomalous field dependence is observed only in 
the clean limitel, suggesting the importance of nonlocal 
effects in understanding the field dependence of j(H). 
Moreover, it has been reported that the vortex core ra- 
dius depends on applied magnetic fiftld and shrinks at 
higher fields in NbSe 2 B and in CeRu 2 Q. 

Another complication especially for borocarbides 
(RNi 2 B 2 C, R = rare earth) is that a square FLL is formed 
in some of these compounds at high magnetic .fields, 
whereas a hexagonal FLL is realized at low fields. ETEI This 
is not expected for the local model with isotropic inter- 



vortex interactions and thereby suggests the importance 
of considering electronic structure (or the Fermi surface) 
and the associated nonlocal corrections in the specific 
compound. 

We report on /xSR measurements of the magnetic field 
dependence of the a-b magnetic penetration depth A, the 
effective vortex core radius p v , and the apex angle of the 
FLL 6 in single crystalline YNi 2 B 2 C. We demonstrate 
that the proper reconstruction of the field profile with a 
square FLL is obtained from the pSK spegtra only when 
the nonlocal corrections are considerred.Q The field de- 
pendence of A turned out to be linear over the entire 
magnetic field range of observation. More importantly, 
it was found that p v shrinks sharply with increasing mag- 
netic field and levels off at higher fields. This shrinkage, 
however, is much steeper than that expected for the case 
when the y/~H behavior of j(H) is entirely attributed to 
that of p v , strongly suggesting that the anomaly in j(H) 
is mostly from the quasiparticle excitations outside the 
vortex cores. 

The single crystal of YNi 2 B 2 C used in this experiment 
(residual resistivity ratio, or rrr ~ 37.4) had a surface 
area of ~ 64 mm 2 . The superconducting transition tem- 
perature T c and the upper critical field H C 2 (T = 3 K) 
determined from resistivity and specific tpat measure- 
ments were 15.4 K and 7.0 T, respectively.! 2 ! ^tSR experi- 
ments were performed on the M15 and M20 surface muon 
beamlines at TRIUMF. An experimental setup with high 
timing resolution was employed to measure the transverse 
field (TF-) /xSR time spectra up to 5 T. The sample was 
mounted with its c-axis parallel to the applied field and 
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beam directions, while the initial muon spin polarization 
was perpendicular to the applied field. The sample was 
field cooled at the measured magnetic fields to minimize 
disorder of the FLL due to flux pinning. Since the muons 
stop randomly on the length scale of the FLL, the muon 
spin precession signal provides a random sampling of the 
internal field distribution in the FLL state. 
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FIG. 1. Fourier Transform of the fiSK time spectra in 
YNi2B2C at 3 K with a strong apodization (see text). 

Figure |l| shows the fast Fourier transforms (FFT) of 
the muon precession signal in YM2B2C fur different fields 
at T ~ 3.0 K with strong apodizationE£l. The real am- 
plitude of the FFT corresponds to the internal magnetic 
field distribution in the FLL state convoluted with an ad- 
ditional damping to account for the weak nuclear dipo- 
lar fields, FLL disorder, and distortions originating from 
the finite time window and the reduced statisties at later 
timesJiil The high-field cutoff reflecting the magnetic field 
at the vortex core is clearly observed. The small peak 
near H — H = is the residual background generated 
by muons which missed the sample. 

In our preliminary analysistHI, it was revealed that the 
local London model with a square FLL fails to reproduce 
the observed /iSR spectra in YNi 2 B 2 C. More specifically, 
the apex angle 9 of the FLL gradually increases from 
60° with increasing field, but it levels off over the field 
range above 0.5 T with 9 ~ 75° (see Fig. |](c)) where the 
square FJLL is-, established by other measurements (i.e., 
9 ~ 90°)LHi3~K Thus, the result in Ref. [l|was obtained 
with 9 ~ 75° for H > 0.5 T. As it is demonstrated below, 
we have found that this problem is. alleviated by taking 
account of the nonlocal corrections!! The local magnetic 
field at any point in the a-b plane is 



H(r) 



H ° 5 1 + K 2 X 2 + A 4 (0.0705Cif 4 + 0.675Cfc2fc2) 

(i) 

where K is the reciprocal lattice vector, 
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FIG. 2. The muon precession signal P x (t) in YN12B2C at 
H = 0.9 T, displayed in a rotating-reference-frame frequancy 
of ~2 MHz. For the solid/dashed curves, see the text. 
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with x and y being the plane of precession, a the FLL 
parameter, 9 the apex angle of the FLL, Hq the average 
magnetic field, A the magnetic penetration depth, and 
being the cutoff parameter. The above reciprocal lat- 
tice vectors correspond to the case where the diagonal 
direction of the FLL (= u + v = 2sin |x) is along the 
(100) direction of the crystal axis. The anisotropic pa- 
rameter C is determined by the band structure, in which 
C scales with A as C = Co/A 2 B The coefficients for C 
in Eq.(|l|) werej-^dopted from the theoretical estimation 
for LUN12B2C.EI The local London model is obtained by 
putting C = 0. 

In addition to the nonlocal corrections, we have de- 
veloped a program to analyze the ^SK spectra in the 
time domain to eliminate the uncertainty in the estima- 
tion of statistical errors associated with fitting the FFT 
spectra. The theoretical time evolution of the muon spin 
polarization was generated by assuming the field profile 
ofEq. 



P x (t)+iPy(t) = 



dr 



dH(r] 



exp(i 7AtJ ff(r)i)cLff (5) 



(7^ is the muon gyromagnetic ratio) and compared with 
the time spectra by the chi-square (x 2 ) minimization 
technique. Considering|-the-,results of small angle neu- 
tron scattering (SANS)0ii3'&3 and scanning tunneling mi- 
croscopy/spectroscopy (STM/STS)B, the apex angle 9 
was fixed to 90° for H > 0.4 T while it was treated as a 
fitting parameter for H < 0.4 T. 

A typical example of the //SR time spectra measured in 
YNI2B2C under a magnetic field of 0.9 T is shown in Fig. 
|2|, where the solid curve is a fit by the nonlocal London 
model while the dashed curve is by the local model with 
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the apex angle fixed to 90°. The value of deduced x 2 
for the nonlocal model is more than three times smaller 
than that for the local model, indicating that the nonlo- 
cal model provides much better description of the data. 
The rate of additional Gaussian relaxation due to trivial 
sources (nuclear dipolar fields, vortex pinnning, etc.) is 
about 0.34 /is -1 at 0.9 T and it tends to be idependent of 
the field. Figure ||| shows the contour plot of H(r) around 
a vortex at H = 0.9 T reproduced from /iSR data, where 
the (100) axis of the crystal is along the horizontal direc- 
tion. The fourfold symmetry due to the nonlocal correc- 
tions in Eq.(0) is clearly observed. We note that there 
are two possible orientations of the FLL configuration in 
Eq.(|l]), where the diagonal direction u + v is parallel with 
either the (100) or the (110) crystalline axis. We have 
found that the field distribution with u + v parallel with 
(110) does not reproduce our data with any combination 
of parameters. This is perfectly in line with the results of 
other experiments, as well as the theoretical calculation 
which^yields a lower free energy for u + v parallel with 
(100)1 
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FIG. 3. The contour map of a flux-line at H — 0.9 T in 
real space, where the unit cell length a of FLL is 471.9 A. 

The physical parameters A, £„, 9 and C versus nor- 
malized external field at 3 K are shown in Fig. ||. 
We treated C as a fitting parameter because its value 
in YNi 2 B 2 C is unknown. The A in YNi 2 B 2 C clearly 
exhibits a linear 7?-dependence. A fit to the relation 
X(h) = A(0)(1 + rj ■ h), (h = H/H c2 ) provides a dimcn- 
sionless parameter rj that represents the strength of the 
pair-breaking effect. We obtain rj = 0.97 (with A(0) = 
567.8 A) which is slightly smaller than that in NbSe 2 (i.e., 
r/ = 1.61 at 0.33TJ). The cutoff parameter (Fig. |(b), 
solid squares) shows a steep decrease with increasing H 
and subsequently levels off at h ^-H/H C 2 > 0.1 (H > 
0.7 T). In our preliminary analysis^, we interpreted this 
cutoff parameter as p v (see eq. (2) in Ref. |l^). In the field 
region h < 0.06 where 9 was set as a free parameter, 9 
gradually decreases with decreasing field, indicating that 
the FLL transforms into a nearly hexagonal lattice. How- 
ever, 9 does not reach 60° in the lowest magnetic field. 
The anisotropy C decreases with increasing H, where it 



exhibits little correlation with 9. While the value at lower 
field is close to the theoretical estimation (~ 0.22 at 0.05 
T in LuNi 2 B 2 CS), we found that CA 2 tends to decrease 
with increasing field. 

The field profile in Fig. || implies that there is an 
anisotropy between the (100) and (110) directions in the 
effective length scales (A and p v ), whereas the model pa- 
rameters in Eq.(Q) represent mean values. Moreover, spe- 
cial precaution must be taken to interpret the parameters 
in Eq.(|l|) upon the introduction of nonlocal corrections 
involving the higher order terms of K, where the defini- 
tion of these length scales are modified from those found 
in the previous analysis with local London models, mak- 
ing it unsuitable to compare directly. In order to evaluate 
p v including the effect of anisotropy, we calculated the 
supercurrent density J(r) from the deduced H(r) using 
Maxwell's relation J(r) = | Vx H(r)|. The radius p v was 
then defined as the distance from the vortex center for 
which J(r) reaches its maximum value. The estimated 
values are pwioo) — 66. 7A and p v /%io) = 70. 8A, yielding 
the ratio p v ( W o)/ Pv(iw) = 0-942. 

The field dependence of p v is shown in Fig. 0(b). The 
values of p v are systematically larger than , suggesting 
that it may not be appropriate to interpret the cutoff 
parameter as the vortex core radius, wiiereas we assumed 
that p v = in the previous analysialj. We stress that 
the core radius can be obtained directly from the field 
profile H(r) deduced from the pSR data, independent of 
the details of the FLL model used.Eil Having said this, 
the field dependence of p v is qualitativly similar to that 
of £„, showing a steep decrease with increasing field in 
the field range H/H c2 < 0.15. 
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FIG. 4. The H dependence of (a) A, (b) p v determined by 
pSK where /9„(ioo) is shown by open squares, p v mo) by solid 
circles and by solid squares, (c) #r-k|y circles and the one 
deduced from the local London modelE3 by triangles, and (d) 
C in the FFL state of YNi 2 B 2 C at 3 K. The dashed curves in 
(b) are described in the text. 
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Recent calculations for s-wave superconductors based 
on the quasiclassical Eilenberger equations predicts a 
shrinkage of p v due to vortex- vortex interactions .113 The 
quasiparticle density of states (DOS) N(H) is propor- 
tional to H$ with (3 = 0.67 at T = in their prediction. 
Provided that all the DOS comes from inside the vortex 
cores, we would expect 

N(H) = N col . c (H)o,7rpl-H^H^ (6) 

where the factor H arises from the number of vortices per 
unit area, and p v oc i?^" 1 )/ 2 . Fitting the field depen- 
dence of p„(ioo) an d Pu(iio) m Fig. ^(b) to the relation 
Pv = poh^-W 2 yields /3 (100> = 0.026, p 0(100> = 24.7A 
and P(no) — 0.014, po(no) = 27. lA, while ±he field de- 
pendence of j(H) yields (3 = 0.430 = (3 SB E. The con- 
siderably smaller values for /3/ioo> and Alio) compared 
with /?sh or the theoretical prediction strongly suggests 
that the origin of the \H behavior of "/(H) is related 
to the quasiparticle excitations outside the vortex cores. 
This is in marked contrast with the case of CeRu2 where 
(3 ~ /3sh, indicating that the DOS is mostly attributed 
to the quasiparticles within the vortex cores.Q The ex- 
istence of delocalized quasiparticle excitations is further 
suggested by the fact that de Haas- van Alphen effect has 
been clearly observed in the mixed state of YM2B2C, 
where the cyclotron radius is much larger than the co- 
herence length £.EZl Surface impedance Z s measurements 
also indicate delocalized quasiparticles outside the vortex 
cores.Ej The magnetic field dependence of N(H) inside 
the cores estimated by Z s is proportional to H, except at 
very low field. These results are consistent with our con- 
clusion that the localized quasiparticles within the vortex 
cores (determined by p v ) contribute little to the \/~H be- 
havior of the Sommerfeld constant, at least for h > 0.15. 
Here, we note that the agreement between (3 and /3sh 
is improved by assuming that N(H) oc p v ■ H instead 
of Eq. (H)t3, although the microscopic origin of this lin- 
ear relation is not obvious at this stage. In any case, 
the small (3 and associated steeper field dependence of p v 
at lower fields might be partly explained by the multi- 
band effect, where the electronic^structure is effectively 
described by a two-band modeled. The RCS coherence 
length £0 = hvp/nAo (where p v < 0.6£o") estimated 
from the Fermi velocity Vp and the energy gap A in 
YNi 2 B 2 C, is 60 to 120 A for one group and 370 A for 
another branch, suggesting that p v (H — > 0) is controlled 
by the larger value of £o- We also note a possible connec- 
tion to the anisotropic energy sap in YM2B2C reported 
by photoemission spectroscopy^, where £0 is scaled by 
the magnitude of Ao. 

Finally, we discuss the apex angle 9 at lower fields 
where the deviation from a hexagonal lattice is expected 
(h < 0.04). We found that the agreement between the 
measured field distribution and calculations based on the 
present model becomes far from satisfactory in the field 
range at ft, < 0.04 (H < 0.3 T). This is probably due to 
the presence of the deep minima along the (110) direction 



in Fig. |3J, which persists irrespective of the apex angle 
(Note that the square shaped field distribution is inde- 
pendent of 9, as is evident in Eq.(|)). The poor agree- 
ment strongly suggests that Eq. (Q) gives the true ground 
state only for the case of a square FLL, while the more 
isotropic distribution would be realized at lower fields,-as 
the 9 ~ 60° well reproduced by the local London modclll3. 
Thus, a more refined model is needed to reproduce the 
complete evolution of the FLL with field. We also point 
out the possibility that FLL domains present through the 
hexagonal-to-square transition (like in LUM2B2O2 2 ]) play 
an important role. 

In summary, we found that p v shrinks steeply with in- 
creasing field while A depends linearly on the magnetic 
field, strongly suggesting the presence of excess quasi- 
particles outside the vortex cores at higher fields. These 
results indicate the need to reconsider the conventional 
picture of a rigid normal-electron core by taking into ac- 
count the vortex-vortex interactions mediated by delo- 
calized quasiparticles. 

We thank the TPJUMF pSR staff for technical sup- 
port. This work was partially supported by a JSPS 
Research Fellowships for Young Scientists, Japan, by a 
Grant-in- Aid for Science Research on Priority Areas from 
the Ministry of Education, Culture, Sports, Science, and 
Technology, Japan and also by a Grant from the CREST, 
JST, Japan. 



' Also at School of Mathematical and Physical Science, The 
Graduate University for Advanced Studies 

1 M. Hedo, et al, J. Phys. Soc. Jpn. 67, 272 (1998). 

2 M. Nohara, et al, J. Phys. Soc. Jpn. 68, 1078 (1999). 

3 J.E. Sonier, et al, Phys. Rev. Lett. 79, 1742 (1997). 

4 R. Kadono, et al, Phys. Rev. B 63, 224520 (2001). 

5 U. Yaron, et al, Nature 382, 236 (1996). 

6 M.R. Eskildsen, et al, Phys. Rev. Lett. 78, 1968 (1997). 

7 D.McK. Paul, et al, Phys. Rev. Lett. 80, 1517 (1998). 

8 M.R. Eskildsen, et al, Nature 393, 242 (1998). 

9 V.G. Kogan, et al, Phys. Rev. B 55, R8693 (1997). 

10 R.H. Norton et al, J. Opt. Sci. Soc. Am. 66, 259 (1976); 
67, 419 (1976). 

11 J.E. Sonier, et al, Rev. Mod. Phys. 72, 769 (2000). 

12 K. Ohishi, et al, Physica B 289-290, 377 (2000). 

13 M. Yethiraj, et al, Phys. Rev. Lett. 78, 4849 (1997). 

14 M. Yethiraj, et al, Phys. Rev. B 58, R14767 (1998). 

15 H. Sakata, et al, Phys. Rev. Lett. 84, 1583 (2000). 

16 M. Ichioka, et al, Phys. Rev. B 59, 184 (1999). 

17 T. Terashima, et al, Phys. Rev. B 56, 5120 (1997). 

18 K. Izawa, et al, Phys. Rev. Lett. 86, 1327 (2001). 

19 J. E. Sonier, et al, Phys. Rev. Lett. 82, 4914 (1999). 

20 S. V. Shulga, et al, Phys. Rev. Lett. 80, 1730 (1998). 

21 T. Yokoya, et al, Phys. Rev. Lett. 85, 4952 (2000). 

22 L.Ya. Vinnikov, et al, Physica B 284-288, 813 (2000) 



4 



